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Abstract: We consider a two-component one-dimensional model of gap 
solitons (GSs), which is based on two nonlinear Schrodinger equations, 
coupled by repulsive XPM (cross-phase-modulation) terms, in the absence 
of the SPM (self-phase-modulation) nonlinearity. The equations include a 
periodic potential acting on both components, thus giving rise to GSs of 
the "symbiotic" type, which exist solely due to the repulsive interaction 
between the two components. The model may be implemented for "holo- 
graphic solitons" in optics, and in binary bosonic or fermionic gases trapped 
in the optical lattice. Fundamental symbiotic GSs are constructed, and their 
stability is investigated, in the first two finite bandgaps of the underlying 
spectrum. Symmetric solitons are destabilized, including their entire family 
in the second bandgap, by symmetry -breaking perturbations above a critical 
value of the total power. Asymmetric solitons of intra-gap and inter-gap 
types are studied too, with the propagation constants of the two components 
falling into the same or different bandgaps, respectively. The increase of the 
asymmetry between the components leads to shrinkage of the stability areas 
of the GSs. Inter-gap GSs are stable only in a strongly asymmetric form, in 
which the first-bandgap component is a dominating one. Intra-gap solitons 
are unstable in the second bandgap. Unstable two-component GSs are 
transformed into persistent breathers. In addition to systematic numerical 
considerations, analytical results are obtained by means of an extended 
("tailed") Thomas-Fermi approximation (TFA). 

©2012 Optical Society of America 

OCIS codes: (190.6135) Spatial solitons; (160.5293) Photonic bandgap materials; (020.1475) 
Bose-Einstein condensates. 
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1. Introduction 

Studies of solitons in spatially periodic (lattice) potentials have grown into a vast area of re- 
search, with profoundly important applications to nonlinear optics, plasmonics, and matter 
waves in quantum gases, as outlined in recent reviews HUD- In ultracold bosonic and fermionic 
gases, periodic potentials can be created, in the form of optical lattices, by coherent laser beams 
illuminating the gas in opposite directions J5H7). Effective lattice potentials for optical waves 
are induced by photonic crystals, which are built as permanent structures by means of vari- 
ous techniques |2, 8]|9), or as laser-induced virtual structures in photorefractive crystals iflOl . 
Parallel to the progress in the experiments, the study of the interplay between the nonlinearity 
and periodic potentials has been an incentive for the rapid developments of theoretical meth- 
ods 01 111121 . Both the experimental and theoretical results reveal that solitons can be created in 
lattice potentials, if they do not exist in the uniform space [this is the case of gap solitons (GSs) 
supported by the self-defocusing nonlinearity, see original works |fl3UT6l and reviews l7l [T7) l, 
and solitons may be stabilized, if they are unstable without the lattice (multidimensional soli- 
tons in the case of self-focusing, as shown in Refs. [18-25), see also reviews OQGDS))- The 
stability of GSs has been studied in detail too — chiefly, close to edges of the corresponding 
bandgaps — in one 11271^291 and two |30l dimensions alike. 

An essential extension of the theme is the study of two-component solitons in lattice po- 
tentials. In particular, if both the self -phase-modulation and cross-phase-modulation (SPM and 
XPM) nonlinearities, i.e., intra- and inter-species interactions, are repulsive, one can construct 
two-component GSs of intra-gap and inter-gap types, with chemical potentials of the compo- 
nents (or propagation constants, in terms of optical media) falling, respectively, into the same 
or different bandgaps of the underlying linear spectrum |[3"Tl[3"2l . In the case of the attractive 
SPM, a family of stable semi-gap solitons was found too, with one component residing in the 
infinite gap, while the other stays in a finite bandgap [32|. The GSs supported by the XPM 
repulsion dominating over the intrinsic (SPM-mediated) attraction may be regarded as an ex- 
ample of symbiotic solitons. In the free space (without the lattice potential), symbiotic solitons 
are supported by the XPM attraction between their two components, despite the action of the 
repulsive SPM in each one [33-35]. This mechanism may be additionally enhanced by the lin- 
ear coupling (interconversion) between the components 136) . Another case of the "symbiosis" 
was reported in Ref. f3~7l , where the action of the lattice potential on a single component was 
sufficient for the stabilization of two-dimensional (2D) two-component solitons against the col- 
lapse, the stabilizing effect of the lattice on the second component being mediated by the XPM 
interaction. In addition, the attraction between the components, competing with the intrinsic 
repulsion, may cause spatial splitting between two components of the GS, as for these compo- 
nents, whose effective masses are negative ifTSI . the attractive interaction potential gives rise to 



a repulsion force 0211381 . 

The ultimate form of the model which gives rise to two-component GSs of the symbiotic 
type is the one with no intra-species nonlinearity, the formation of the GSs being accounted 
for by the interplay of the repulsion between the components and the lattice potential acting 
on both of them. In optics, the setting with the XPM-only interactions is known in the form of 
the "holographic nonlinearity", which can be induced in photorefractive crystals for a pair of 
coherent beams with a small angle between their wave vectors, giving rise to single- 13911401 
and double-peak fiTl solitons. Both beams are made by splitting a single laser signal, hence 
the power ratio between them (which is essential for the analysis reported below) can be var- 
ied by changing the splitting conditions. The creation of 2D spatial "holographic solitons" in 
a photorefractive-photovoltaic crystal with the self-focusing nonlinearity was demonstrated in 
Ref. l42ll (such solitons are stable, as the collapse is arrested by the saturation of the self- 
focusing). To implement the situation considered here, the sign of the nonlinearity may be 
switched to self-defocusing by the reversal of the bias voltage, and the effective lattice po- 
tential may be induced by implanting appropriate dopants, with the concentration periodically 
modulated in one direction, which will render the setting quasi-one-dimensional. 

In binary bosonic gases, a similar setting may be realized by switching off the SPM nonlin- 
earity with the help of the Feshbach resonance, although one may need to apply two different 
spatially uniform control fields (one magnetic and one optical) to do it simultaneously in both 
components. On the other hand, the same setting is natural for a mixture of two fermionic com- 
ponents with the repulsive interaction between them, which may represent two states of the 
same atomic species, with different values of the total atomic spin {F). If spins of both compo- 
nents are polarized by an external magnetic field, the SPM nonlinearity will be completely sup- 
pressed by the Pauli blockade while the inter-component interaction remains active [6], hence 
the setting may be described by a pair of Schrodinger equations for the two wave functions, 
coupled by XPM terms. 

The objective of this work is to present basic families of one-dimensional symbiotic GSs, 
supported solely by the repulsive XPM nonlinearity in the combination with the lattice poten- 
tial, and analyze their stability, via the computation of eigenvalues for small perturbations and 
direct simulations of the perturbed evolution. The difference from the previously studied mod- 
els of symbiotic solitons 1311321 is that the solitons where created there by the SPM nonlinearity 
separately in each component, while the XPM interaction determined the interaction between 
them and a possibility of creating two-component bound states. Here, the two-component GSs 
may exist solely due to the repulsive XPM interactions between the components. 

We conclude that the symmetric solitons, built of equal components, are destabilized by 
symmetry-breaking perturbations above a certain critical value of the soliton's power. The 
analysis is chiefly focused on asymmetric symbiotic GSs, and on breathers into which un- 
stable solitons are transformed. The model is introduced in Section II, which is followed by 
the analytical approximation presented in Section III. It is an extended ("tailed") version of 
the Thomas-Fermi approximation, TFA, which may be applied to other models too. In Section 
IV, we report systematic numerical results obtained for fundamental solitons of both the intra- 
and inter-gap types, hosted by the first two finite bandgaps of the system's spectrum. The most 
essential findings are summarized in the form of plots showing the change of the GS stability 
region with the variation of the degree of asymmetry of the two-component symbiotic solitons, 
which is a new feature exhibited by the present system. In particular, the stability area of intra- 
gap solitons shrinks with the increase of the asymmetry, while inter-gap solitons may be stable 
only if the asymmetry is large enough, in favor of the first-bandgap component, and intra-gap 
solitons in the second bandgap are completely unstable. The paper is concluded by Section V. 
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Fig. 1. The identical bandgap structures produced by the linearization of Eqs. © and (O 
for e = 6. Shaded areas are occupied by the Bloch bands, where gap solitons do not exist. 



2. The model 

The model outlined above is represented by the system of XPM-coupled Schrodinger equations 
for local amplitudes of co-propagating electromagnetic waves in the planar optical waveguide, 
u(x,z) and v(x,z), where x and z are the transverse coordinate and propagation distance, without 
the SPM terms, and with the lattice potential of depth 2e > acting on both components: 

du 1 d 2 u , l9 , . 

«^- + ^^rr- \v M + ecos (2x )u = 0, (1) 
dz 2 dx L 

i^ + ^^rZ-\u\ 2 v + ecos(2x)v = 0. (2) 
dz 2 dx z 

The variables are scaled so as to make the lattice period equal to n, and the coefficients in front 
of the diffraction and XPM terms equal to 1. In the case of matter waves, u and v are wave 
functions of the two components, and z is replaced by time t . Direct simulations of Eqs. (Q]i and 
(O were performed with the help of the split-step Fourier-transform technique. 
Stationary solutions to Eqs. ([TJ, <J3J are looked as 

u(x,z) =e ikz U(x) 1 v(x,z)=e iqz V(x), (3) 

where the real propagation constants, k and q, are different, in the general case, and real func- 
tions U (x) and V (x) obey equations 

-kU + hj" -V 2 U + ecos(2x)U = 0, (4) 

-qV + -V" - U 2 V + e cos{2x)V = 0, (5) 

with the prime standing for djdx. Numerical solutions to Eqs. (0]i and © were obtained by 
means of the Newton's method. 

Solitons are characterized by the total power, 

\\U\ 2 + \V\ 2 )dx = P u + P v , (6) 



with both P„ and P v being dynamical invariants of Eqs. ([T), @, and by the asymmetry ratio, 



R = (P U -P V )/(P U + P V ). (7) 

The total power and asymmetry may be naturally considered as functions of the propagation 
constants, k and q. 

The well-known bandgap spectrum of the linearized version of Eqs. @, © (see, e.g., book 
ifTTl ) is displayed in Fig. [TJ the right edge of the first finite bandgap being 

= 6) » 3.75. (8) 

The location of GSs is identified with respect to bandgaps of the spectrum. In this work, results 
are reported for composite GSs whose two components belong to the first and second finite 
bandgaps. 

Stability of the stationary solutions can be investigated by means of the linearization against 
small perturbations |26-28|. To this end, perturbed solutions of Eqs. ([]]) and (O are looked for 

as 



u(x,z) = e ikz U(x) + u 1 (x)e- az + u* 2 (x)e iX * 



v(x,z) = e"< z V{x)+ Vl (x)e- IAZ + v* 2 (x) 



(9) 



where u\ t2 and vi 2 are wave functions of infinitesimal perturbations, and X is the respective 
instability growth rate, which may be complex (the asterisk stands for the complex conjugate). 
The instability takes place if there is at least one eigenvalue with Im(A) > 0. The substitution 
of ansatz Eqs. (0 into Eqs. £[), (f2j) and the linearization with respect to the small perturbations 
leads to the eigenvalue problem based on the following equations: 



av i - ~v" + U 2 (x)vi+U(x)V(x)(ui+u 2 )-ecos(2x)vi 





(10) 


Av 2 , 


(ID 


Ami, 


(12) 


Xu 2 - 


(13) 



-qv 2 + ivjj - U 2 (x)v 2 -U(x)V(x)(u 1 + u 2 ) + £cos(2jc)v 2 

ku\ — —u'{ + V 2 (x)u-{ +U(x)V(x) (vi +V2) — ecos(2x)t/i 

-ku 2 + ^u' 2 ' - V 2 (x)u 2 - U(x)V(x) (vi + v 2 ) + £cos(2x)m 2 

These equations were solved by means of the fourth-order center-difference numerical scheme. 

Results for the shape and stability of GSs of different types are presented below for lattice 
strength e = 6, which adequately represents the generic case. Note, in particular, that Fig. [T] was 
plotted for this value of the lattice-potential's strength. 

3. The extended Thomas-Fermi approximation 

It is well known that, close to edges of the bandgap, GSs feature an undulating shape, which 
may be approximated by a Bloch wave function modulated by a slowly varying envelope 
lfl3l[T5l . On the other hand, deeper inside the bandgap, the GSs are strongly localized (see, 
e.g., Figs. H] and [7] below), which suggests to approximate them by means of the variational 
method based on the Gaussian ansatz 13 21 1431 . This approximation was quite efficient for the 
description of GSs in single-component models l43l . while for two-component systems it be- 
comes cumbersome 113 1 11321 . 

Explicit analytical results for well-localized patterns can be obtained by means of the TFA 
lO, which, in the simplest case, neglects the kinetic-energy terms, U" and V" ', in Eqs. ©, 
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Fig. 2. (a) The continuous (blue) curves show the numerically found amplitude of the fun- 
damental symmetric gap solitons (with equal components), versus propagation constant 
k = q, in the first and second bandgaps at e = 6.0. The chain of symbols is the analytical 
approximation for the same dependence, as produced by the improved TFA in the form 
of Eq. 1161) . The dashed curve is the result of the usual TFA, which corresponds to Eq. 
i Kit without the correction (second) terms, (b) Total power P for the same soliton fami- 
lies, whose stable and unstable portions are designated by the bold dotted and dashed lines, 
respectively. The latter one is destabilized by symmetry-breaking perturbations, while the 
entire family is stable in the framework of the single-component model. Coordinates of the 
stability/instability border are given by Eq. J2U . The chain of squares shows the analytical 
dependence produced by the TFA, see Eq. d20t . 



Assuming, for the sake of the definiteness, q < k and also \k\ < e (the TFA is irrelevant for 
\k\ > e), the approximation yields the fields inside the inner layer of the solution: 



U 2 (x) 
V 2 (x) 



e cos (2x) — q 
e cos (2x) — k 



at \x\ < xq : 



(14) 



Thus, the TFA predicts the core part of the solution in the form of peaks in the two components 
with the same width, 2xq, but different heights, {U 2 7 V 2 } max = {E — q,E — k}. This structure 
complies with numerically generated examples of asymmetric solitons displayed below in Figs. 

Ekb). 

Further, the expansion of expressions ( [Pfl ) around the soliton's center (x = 0), yields 



U{x) 
V(x) 



y /£-g - (e/ ye-gj s 2 
y/e — k— (e/\/e-k)x 2 



(15) 



The substitution of the second derivatives of the fields at x = 0, calculated as per Eq. ( fl5l ), into 
Eqs. dUi, (O yields a corrected expression for the soliton's amplitudes: 



U(x = 0) 
V(x = 0) 



x/e^q - e [2 (e - k) yfe=q\ 



-l 



[ y/e-k-£[2{£-q)y/e-k\~ 
along with the condition for the applicability of the TFA: 

e < (e - k) (e-q). 



(16) 



(17) 



For the symmetric-GS families in the two first finite bandgaps, the amplitude predicted by the 
improved TFA in the form of Eq. dT6b is displayed, as a function of k = q, in Fig. |2a) and 
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Fig. 3. The evolution of a weakly unstable single-component fundamental soliton at k = 
-3.7. 



compared to its numerically found counterpart. It is worthy to note that the correction terms in 
Eq. ( [Tol l essentially improve the agreement of the TFA prediction with the numerical findings. 

At |jc| > Xo, the TFA gives V(x) = 0, which is a continuous extension of the respective ex- 
pression Eq. O in the V component, while the continuity of fields U (x) and U'(x) makes it 
necessary to match the respective expression Eq. ( TPfl i to "tails", which, in the lowest approxi- 
mation, satisfy equation U" = 0. The continuity is provided by the following tail solution: 



Vk~q- yj(e 2 -k 2 )/(k-q)(\x\ -Xq) 
at < |jc| — jcd <(k-q) /Ve 2 -k 2 ; 
0, at \x\ -xo>(k-q) /Ve 2 -k 2 . 



(18) 



The integration of expressions Eq. (fT4] i and Eq. (fT8l yields the following approximation for 
the powers of the two components: 



Pu 



TFA 



Ve 2 — k 2 -qcos l {k/e) 
+ (2/3) {k -qf/y/e^l?; 
Ve 2 -k 2 -kcos- l (k/e). 



(19) 



The substitution of approximation Eqs. $1% into definitions Eq. © and Eq. ^ of the total 
power and asymmetry demonstrates an agreement with numerical results. For instance, the 
slope of the curve R(q) for the intra-gap GSs at fixed k (see Fig.[8]below) at the symmetry point 
(k = q), as predicted by Eqs. ( fT9b for e = 6 and k = 1, is (dR/dq) \ q= i l m —0.155, while its nu- 
merically found counterpart is w —0. 160. Further, the analysis of Eqs. (fT9b readily demonstrates 
that the strongly asymmetric solitons may exist up to the limit ofR—> 1, which is corroborated 
by the existence area for the intra-gap solitons shown below in Fig.rTTTb). 

Another corollary of Eqs. ( fT9l is the prediction for the total power for the symmetric solitons, 



P(k = q)=2 y / e 2 -k 2 -kcos- 1 (k/e) 



(20) 



which is plotted in Fig.|2jb), along with its numerically found counterpart. Although the TFA 
does not predict edges of the bandgaps, the overall analytical prediction for P(k) runs quite close 
to the numerical curve, except for near the right edge, where, indeed condition Eq. ( fTTT i does 
not hold for e = 6 and k = 3.75, see Eq. ©. Note that this very simple analytical approximation 
was not derived before in numerous works dealing with single-component GSs. 
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Fig. 4. Examples of fundamental symmetric gap solitons found in the first and second 
bandgaps, for k = q = 2.0 and k = q = —2.0 (left and right panels, respectively). Here and 
in similar figures below, the background pattern (green sinusoid) represents the underly- 
ing periodic potential. Both solitons are stable as solutions of the single-component model, 
but only the one corresponding to k = q = 2.0 remains stable in the two-component sys- 
tem, while its counterpart pertaining to k = q = — 2 is destabilized by symmetry-breaking 
perturbations. 
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Fig. 5. (a) and (b) The spontaneous transformation of an unstable symmetric fundamental 
soliton for u- and v-component, in the first bandgap, with k = q = 0, into a stable asym- 
metric breather, (c) The top and bottom plots display, respectively, the evolution of the 
peak-power difference, max(|u (jc,z)| 2 ) — max(\u(x,z)\ 2 ), and the separation between cen- 
ters of the two components, X u and X v , which as per Eq. J22b . 
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Fig. 7. Examples of stable solitons of the intra-gap type found in the first finite bandgap, 
with fixed asymmetry R = -0.5: (a) k = 3 and q = 3.4601; (b) k = 1 and q = 2.843; (c) 
k = —0.5 and q = 2.51 16. Fields U (x) and V(x), which pertain to propagation constants k 
and q, are shown, respectively, by the magenta (lower) and blue (higher) profiles. 



4. Results of the numerical analysis 



4.1. Symmetric solitons 

Obviously, the shape of symmetric solitons, built of two equal components [with k = q and 
U(x) = V(x), see Eq. (0)], is identical to that of GSs in the single-component model. How- 
ever, there is a drastic difference in the stability of the symmetric GSs between the single- 
and two-component systems. Almost the entire symmetric family is stable against symmetric 
perturbations, i.e., it is stable in the framework of the single-component equation (in agree- 
ment with previously known results ifTTI '). except for a weak oscillatory instability, accounted 
for by quartets of complex-conjugate eigenvalues, in the form of A = ±/Im (A) ± Re (A) (with 
two mutually independent signs ±), which appears near the left edge of the second bandgap — 
namely, at k < £ m ; n w —3.45. An example the development of the latter instability is displayed 
below in Fig. [3] 

On the other hand, Fig.|2]demonstrates that a considerable part of the family in the first finite 
bandgap, and the entire family in the second bandgap are unstable against symmetry-breaking 
perturbations in the two-component system. The boundary separating the stable and unstable 
subfamilies of the fundamental symmetric GSs in the first finite bandgap corresponds to the 
power and propagation constants is found at 

Per W 7.19, Act « 1.05, (21) 

the symmetric solitons being stable in the intervals of < P < 7 . 1 9, 1 .05 < A: < £ max ~ 3 .75 [see 
Eq. ©]. These results were produced by a numerical solution of Eqs. (TTOb - (TT3T> (the instability 
is oscillatory, characterized by complex eigenvalues). 

Typical examples of stable and unstable fundamental symmetric GSs, found in the first and 
second bandgaps (not too close to their edges), are displayed in Fig. [4] Further, direct simula- 
tions demonstrate that the evolution transforms the unstable symmetric solitons into persistent 
localized breathers, as shown in Figs.|5]and|6] in accordance with the fact that the corresponding 
instability eigenvalues are complex. Although the emerging breather keeps the value of R = 0, 
see Eq. ©, the u- and v- components of the breather generated by the symmetry-breaking 
instability are no longer mutually identical. This manifestation of the symmetry-breaking insta- 
bility is illustrated by Fig.|5jc), which displays the evolution of the difference between the peak 
powers of the two components, and the separation between their centers. The latter is defined 
as 

X u -X v = — f \u(x,z)\ 2 xdx~ — [ \v(x,z)\ 2 xdx. (22) 

Pit J — CO P v J — CO 

It is relevant to mention that the second finite bandgap also contains a branch of the so-called 
subfundamental solitons, whose power is smaller than that of the fundamental GSs 144H46I . 
These are odd modes, squeezed, essentially, into a single cell of the underlying lattice potential. 
The subfundamental solitons are unstable, tending to rearrange themselves into fundamental 
ones belonging to the first finite bandgap, therefore they are not considered below. 

4.2. Asymmetric solitons of the intra-gap type 

As said above, two-component asymmetric fundamental GSs, with different propagation con- 
stants, k ^ q, may be naturally classified as solitons of the intra- and inter-gap types if k and 
q belong to the same or different finite bandgaps 1311 . In this subsection, we report results for 
asymmetric intra-gap solitons with both k and q falling into the first finite bandgap, as well as 
for asymmetric breathers developing from such solitons when they are unstable. 

Examples of stable asymmetric GSs of the intra-gap type are displayed in Fig. [71 for a fixed 
asymmetry ratio, R = —0.5, defined as per Eq. (0. The GS family, along with the family of 
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Fig. 8. The asymmetry ratio, R [defined as per Eq. (7)], versus propagation constant q, at 
fixed values of k = 3.0, 1.0, and k = —0.5 (the top, middle, and bottom curves, respectively), 
for asymmetric fundamental solitons of the intra-gap type. Stable and unstable branches are 
shown by solid and dashed lines, respectively. 
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Fig. 9. An example of a stable strongly asymmetric soliton with k = —0.5 and q = 3.65. 
Fields U (x) and V(x), which pertain to propagation constants k and q, are shown, respec- 
tively, by the magenta (lower) and blue (taller) profiles. 




Fig. 10. A typical example of the transformation of the unstable asymmetric gap solitons 
into a breather, for k = —0.5 and q = 2.0. 
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Fig. 11. (a) The stability border in the plane of the propagation constants, (k,q), for asym- 
metric solitons of the intra-gap type. Only half of the plane is shown, delineated by the 
dotted triangle, within which wavenumbers k and q belong to the first finite bandgap, as 
the other half is a mirror image of the displayed one. (b) The same in the plane of the total 
power and asymmetry ratio, {P,R), defined as per Eqs. l[6) and Q. Localized modes do not 
exist above the right boundary of the stability regions in panel (b). The diagram at R < is 
a mirror image of the one displayed here for R > 0. 



persistent breathers into which unstable solitons are spontaneously transformed, is represented 
in Fig. [8]by dependences R(q) at different fixed values of the other propagation constant, k. 

It is possible to explain the fact that all the R(q) curves converge to R = — 1, as q approaches 
the right edge of the bandgap in Fig. [8] In this case, the V component turns into the delocalized 
Bloch wave function with a diverging power, P v , that corresponds to P u /P v -4 [it is tantamount 
to R — >• — 1, as per Eq. ©]■ An example of a stable GS, close to this limit, with k = —0.5, 
q = 3.65 and/? = —0.981 1, is shown in Fig. [9] The central core of the V-componentis described 
by the TFA, based on Eq. (fT4l i, as the corresponding necessary condition ( fTTl ) holds in this case, 
while the TFA does not apply to the U -component. The presence of undulating tails, which are 
close to the Bloch functions, rather than the simple approximation d 1 St . which is valid far from 
the edge of the bandgap, is also visible in Fig. [8] 

Those asymmetric intra-gap GSs which form unstable subfamilies in Fig. [8] are destabilized 
by oscillatory perturbations. The instability transform the solitons into breathers, see a typical 
example in Fig.[l0][cf. the examples of the destabilization of the symmetric GSs shown in Fig. 
I2b,c)]. The emerging breathers keep values of the asymmetry ratio (j7]i almost identical to those 
of their parent GSs; for instance, in the case displayed in this figure, the unstable soliton with 
^initial = —0.3617 evolves into the breather with Rfi na \ = —0.3623. 

It is relevant to stress that the transformation of unstable stationary GSs into the breathers 
gives rise to little radiation loss of the total power, P. On the other hand, in the general case a 
given unstable gap soliton does not have a stable counterpart with a close value of P, hence this 
unstable soliton cannot transform itself into a slightly excited state of another stable GS. Thus, 
the breathers represent a distinct species of localized modes. 

The most essential results of the stability analysis for the asymmetric solitons of the intra-gap 
type, and for breathers replacing unstable solitons, are summarized by diagrams in the planes 
of {k,q) and (P,R), which are displayed in Fig.[TT] The predictions of the analysis based on 
the computation of the stability eigenvalues for the stationary solitons, as per Eqs. (fT0b-([T3b. 
always comply with stability tests provided by direct simulations of Eqs. (Q]) and (O. 
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Fig. 12. An example of a stable inter-gap soliton, for k = 3 and q = — 1 .65. The single-peak 
and split-peak profiles, (J (x) and V(x), represent, respectively, the components in the first 
and second finite bandgaps. 

As mentioned above, the instability of a part of the branch of the symmetric solitons along 
the line of R = in Fig. fTTT b) implies that the symmetry -breaking perturbations destabilize the 
symmetric solitons in the first finite bandgap at P > P cr , see Eq. d2Tl i, while their counterparts 
are stable in the single-component system. Another clear conclusion is that the stability region 
gradually shrinks with the increase of the asymmetry. 

4.3. Solitons of the inter-gap type 

All the GSs of the inter-gap type, with two propagation constants belonging to the two different 
finite bandgaps, are naturally asymmetric, even if their components have equal powers. Ex- 
amples of stable and unstable inter-gap solitons are displayed in Figs.[T2]and[T3] respectively. 
A noteworthy feature exhibited by these examples is a split-peak structure of the component 
belonging to the second finite bandgap. 

The asymmetry measure, R(q), for families of the inter-gap solitons is plotted in Fig. fPfl' a) 
versus the propagation constant q in the second bandgap, at fixed values of k (the propagation 
constant in the first bandgap). The stability of the respective GS families is also shown in Fig. 
El 

The (in)stability of the inter-gap solitons is summarized by the diagrams in the planes of 
(k,q) and (P,R) presented in Fig. [15] cf. similar diagrams for intra-gap solitons shown above in 
Fig-El As well as in that case, unstable inter-gap solitons are spontaneously replaced by robust 
localized breathers. The spontaneous transformation increases the initial degree of the asym- 
metry: For instance, an unstable inter-gap soliton with R = 0.019 is converted into a breather 
with R = 0.028. 

In the present case too, the existence region of stable modes shrinks with the increase of the 
asymmetry; note also that the stationary inter-gap solitons may be stable solely at sufficiently 
large values of the asymmetry, R > R m i n ss 0.5. The asymmetric shape of the stability diagram 
in Fig.fTBTb) with respect to R > and R < [unlike the symmetry of the diagram for the intra- 
gap solitons implied in Fig. fTlT b)! is explained by the fact that, in definition (O, R > implies 
that the dominant component resides in the first finite bandgap, where it is more robust than in 
the second bandgap. 

Finally, out additional analysis has demonstrated that all the stationary GSs — not only the 
symmetric ones [see Fig. [2), but also all the asymmetric solitons of the intra-gap type — are 
completely unstable in the second finite bandgap. They too tend to spontaneously rearrange 
themselves into breathers, which is not shown here in detail. 




(c) 



Fig. 13. A typical example of the transformation of an unstable intergap soliton into a stable 
breather, at k = and q = —2. (a,b) The evolution of |m| 2 and |v| 2 . (c) The initial profiles 
of U (x) and V(x) (single-peak and split-peak shapes, respectively). 
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Fig. 14. Asymmetry ratio R for the inter-gap solitons versus propagation constant q in the 
second finite bandgap, at fixed values k = 3.0, 2.5, 1.0, and —0.5 (from the top to the 
bottom) of wavenumber k in the first bandgap. Solid and dashed lines designate stable 
stationary solitons and breathers, respectively. 
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Fig. 15. The same as in Fig. QT| but for inter-gap solitons. In (a), the dotted rectangle 
delineates the region occupied by wavenumbers k and q belonging to the first and second 
finite bandgaps, respectively. 



5. Conclusion 

We have introduced the model of symbiotic two-component GSs (gap solitons), based on two 
nonlinear Schrodinger equations coupled by the repulsive XPM terms and including the lat- 
tice potential acting on both components, in the absence of the SPM nonlinearity. The model 
has a realization in optics, in terms of "holographic solitons" in photonic crystals, and as a 
model of binary quantum gases (in particular, a fully polarized fermionic one) loaded into the 
optical-lattice potential. Families of fundamental asymmetric GSs have been constructed in the 
two lowest finite bandgaps, including the modes of both the intra-gap and inter-gap types, i.e., 
those with the propagation constants of the two components belonging to the same or different 
bandgaps, respectively. The existence and stability regions of the symbiotic GSs and breathers, 
into which unstable solitons are transformed, have been identified. A noteworthy finding is that 
symmetry-breaking perturbations destabilize the symmetric GSs in the first finite bandgap, if 
their total power exceeds the critical value given by Eq. (ETT i. along with all the symmetric 
solitons in the second bandgaps. It was demonstrated too that the stability area for the intra- 
gap GSs shrinks with the increase of the asymmetry ratio, R. On the other hand, inter-gap GSs 
may be stable only for sufficiently large ratio, R > 0.5. The intra-gap solitons are completely 
unstable in the second bandgap. Some features of the GS families were explained by means of 
the extended TEA (Thomas-Fermi approximation), augmented by the tails attached to the taller 
component, in the case of asymmetric solitons. 

A natural extension of the analysis may deal with 2D symbiotic gap solitons, supported by the 
square- or radial-lattice potentials. In that case, it may be interesting to consider two-component 
solitary vortices too. 
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